Multisite Entanglement acts as a Better Indicator of Quantum Phase Transitions 

in Spin Models with Three-spin Interactions 
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We employ a genuine multipartite entanglement measure, the generalized geometric measure, for 
investigating the quantum phase transition in an infinite quantum spin- 1/2 chain with two-spin as 
well as three-spin interactions. We show that in contrast to bipartite entanglement of the ground 
state, the genuine multiparty one unerringly indicates the quantum phase transition in the system. 
The system lends itself to a complementary behavior between bipartite and genuine multipartite 
entanglements. We also investigate the trend of the multipartite entanglement in finite chains, in 
which we perform a finite-size scaling to confirm the signature of the quantum phase transition 
already in finite systems. The finite-size investigation renders the phenomenon to be observable in 
real systems. 



I. INTRODUCTION 

Quantum entanglement is one of the basic ingredients 
in quantum physics and is identified as the key resource 
for many potential applications in quantum communica- 
tion and computational tasks [1]. In recent years, great 
efforts have been put towards characterizing and quan- 
tifying the entanglement properties in quantum systems. 
However, a relatively good understanding is limited to 
bipartite systems, with a general characterization and 
quantification of entanglement in multiparty systems be- 
ing still an open problem. 

On the other hand, it has been established in the last 
decade or so that there are connections between bipar- 
tite entanglement [1] and quantum critical phenomena 
[2] in many-body systems. In particular, it has been 
shown that bipartite entanglement can be useful for de- 
tecting (zero temperature) quantum phase transitions, 
which are transitions driven solely by quantum fluctua- 
tions between qualitatively distinct phases of quantum 
many-body systems [3, 4]. The bipartite entanglement 
measures, employed for this purpose, include concurrence 
[5] and logarithmic negativity [6]. 

With the vast developments in the field of simulat- 
ing quantum many-body systems in different physical 
substrates, including atoms in an optical lattice [3, 7], 
trapped ions [8], photons [9], and nuclear magnetic res- 
onance [10], it is now possible to manipulate and mea- 
sure single party, bipartite, and multipartite properties 
of many-body quantum systems in the laboratory. This 
is performed by using the unprecedented levels of con- 
trol that has been recently acquired in preparation and 
manipulation of such systems. Most of these investiga- 
tions are restricted to Hamiltonians with two-body inter- 
actions. 

Quantum many-body systems, however possess, in 
general, also three- and higher-body interactions, along 
with the two-body ones. Three-body interactions give 
rise to a whole new world of interesting physics [11, 12]. It 
has been shown that the multi-body interaction can lead 
to novel quantum phases with striking features [13]. Such 
many-body interactions have been observed, for example, 



in complex oxides [14], as well as in engineered physical 
systems in the laboratory [15-18]. In particular, it was 
proposed that dominant three-body interactions can be 
created in a system of cold polar molecules [16]. Genera- 
tion of effective three-body interactions was also shown to 
be possible by using the two-body collisions of atoms of a 
three-dimensional optical lattice [17]. Three-body inter- 
actions can potentially be simulated in one-dimensional 
atomic lattice gases with Rydberg atoms [18]. Recently, 
Hamiltonians with multi-body interactions in ultracold 
atoms have been realized [19]. 

In this paper, we consider the infinite spin-1/2 chain 
with three-body interactions along with isotropic (two- 
body) XY couplings. The entanglement properties of 
this model have been investigated in Rcfs. [20, 21]. In 
particular, it has been observed that the (bipartite) en- 
tanglement of nearest-neighbor spins cannot reliably in- 
dicate the critical behavior of this model [20] . We show 
that a genuine multiparty entanglement measure of the 
ground state can be used to faithfully indicate the quan- 
tum phase transition (QPT) present. The genuine mul- 
tiparty entanglement employed is the generalized geo- 
metric measure (GGM) [22]. The calculations for the 
infinite chain are carried out analytically, by applying 
the Jordan- Wigncr transformation on the spin Hamilto- 
nian. Keeping in mind that an experimental verification 
of the presented results and effects can potentially be 
carried out in finite spin systems, we have investigated 
the status of the effects, observed for the infinite chain, 
in finite-sized chains also. A finite-size scaling analysis 
reveals that the effect observed in the infinite chain can 
also be reproduced in finite chains. 

The paper is organized as follows. We first briefly 
introduce the genuine multipartite entanglement, the 
GGM, in Sec. II. In Sec. Ill, we present the Hamilto- 
nian of the spin-1/2 XX (i.e., isotropic XY) chain with 
three-spin interactions, that we are investigating. Sec. 
IV presents the results obtained for the infinite and fi- 
nite spin chains. Finally, we give a conclusion in Sec. 
V. 
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II. 



GENERALIZED GEOMETRIC MEASURE 



In this section, we introduce the genuine multiparticle 
entanglement measure which we shall be using to inves- 
tigate the multipartite entanglement of the system under 
consideration. A multiparty pure quantum state is said 
to be genuinely multiparty entangled if it is entangled 
across every bipartition of its constituent parties. To 
quantify genuine multipartite entanglement present in a 
multiparty system, we will use the generalized geometric 
measure [22]. For an N- party pure quantum state \iPn), 
the GGM is defined as 



G{\^ N )) = l-max|(0 JV |V'iv} | s 



(1) 



rr z — lA r - 1 
whence the Hamiltonian takes the form 



JV 



\ (c«c„+i + h.c.) + (c^Cn+2 - h.c.) 



(4) 

The boundary condition employed here is either periodic 
or anti-periodic, depending on whether the number of 
spinlcss fermions is even or odd, respectively. We now 
apply the Fourier transformation 



C n = H ° XP HM4- C n = -j= CX P 



(ikn)ck 



where the maximization is carried out over all pure states 
\4>n) that are not genuinely N-party entangled. The 
above expression can be simplified as 



G(|Vjv» = l-max{A 2 }, 



(2) 



where {A 2 } denotes the collection of the squares of all 
maximal Schmidt coefficients of all bipartite splits of the 
pure state | Vjv) • 



where for an odd (even) number of spinlcss fermions, k = 
(2ir/N)m (k = (2n/N)(m + 1/2)) in the subspaces with 
to = -(JV - l)/2, — (JV - l)/2 + 1, (JV - l)/2 (to = 
-N/2, -N/2 + 1, ...,N/2 - 1). As a consequence, the 
Hamiltonian reduces to 



H 



k 



e(k)c\c k . 



(5) 



The resulting energy dispersion is given by 



III. MODEL HAMILTONIAN 

The model Hamiltonian that we consider in this paper 
has two-body as well as three-body interactions. The 
two-body part is the well-known XX interaction, and is 
interesting on its own [2] . With the introduction of the 
additional three-body terms, the model exhibits interest- 
ing physics, including chiral phases [11]. The governing 
Hamiltonian is given by [23, 24] 



H = - 



J 



N 

[«+l+ a >n+V 



n=l 



(3) 



a 

2 U 



n-l u n°n+l 



where J is the coupling constant for the two-body inter- 
actions while aJ/2 is that for the three-body ones, N is 
the total number of spins arranged in an one-dimensional 
array, and a% (a = x,y,z) are the Pauli spin matrices 
at the site n. Note that a is a dimcnsionlcss quantity 
characterizing the strength of the three-spin interaction. 
Here the periodic boundary condition is assumed i.e., 
cjv+i = fi- This Hamiltonian can be diagonalized in 
two steps. To that aim, we first apply a Jordan- Wigner 
transformation [25] given by 



n 



1 n— 1 

n 



2c m c m j 



2 c m c m 



e(k) = — J cos (k) — — sin (2k) 



(6) 



The e(k) has the following features: In the regime where 
a ^ ol c = 1, only two Fermi points at kp — ±n/2 
are present. But in the regime where a > a c , there 
emerges two additional Fermi points in the negative- 
energy region in the k-space. So, four Fermi points ap- 
pear with the two additional points at kp = arcsin(l/a) 
and kp = it — arcsin(l/o;). In the thermodynamic limit 
(JV —> oo), the ground state of the system corresponds to 
the configuration where all the states with e(k) ^ are 
populated and the s(k) > are empty. So, when a ^ a c , 
the configuration of the ground state is of the states with 
|fc| < kp. On the other hand, when a > a c , the ground 
state is formed by the states with — kp < k < kp and 
kp < k < kp as filled. The energy per site of the ground 
state is 



E o/N = -j- / 
27T J{k} 



e(k)dk, 



(7) 



where the integration region {k} = [— kp,kp] for a ^ 1 
and {k} = [—kp, kp] U [kp, kp] for a > 1. After a simple 
calculation, it takes the form 



E Q /JN = 



1 

7T 
1 

2^ 



a < 1, 



a > 1. 



The ground state energy has a constant value for a < 1 
while it increases with a for a > 1, and the QPT occurs 
at a — a c — 1 [23]. The magnetization for the ground 
state is zero for all values of a. 
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To evaluate the entanglement properties, one needs to 
fully characterize the ground state of the system. Once 
the expectation values of the two-point correlation func- 
tions (cJjC m ) = fn.m are known, the ground state of the 
model can be completely characterized, and any other 
expectation value can be evaluated by using the Wick's 
theorem [25]. The / n)7rl 's for an infinite chain can be 
calculated to form a correlation matrix F, which charac- 
terizes the density matrix of the ground state [21]: 



fn 



(c n C m ) 



— I cxp [— ik(m — n)] dk, (8) 
2tt 7{fc} 



where the integration is again carried out in the region 
{k} = [~k F ,k F ] for a^l and {k} = {-k F ,k F ](j[k F ,k F } 
for a > 1. The analytical expression of the two-point 
correlation function is given by 



1 



/ 77171" \ 



fn,' 



sm 1 

7777T \ 2 / 



1 



2m7r 
for m 7^ and 



[1- (-1H sin 



fn.n 2 " 



to arcsm | — 
a 



a < 1, 

a ^ 1, 
(9) 

(10) 



When 777 is an even integer, the f n ,n+m always vanishes. 
The elimination of the rows and columns from the F- 
matrix corresponds to the spins which do not belong to 
the block, and hence the 77-particle reduced density ma- 
trices (which we denote as p n ) of the ground state can be 
characterized by using the reduced correlation matrix 



F n = 



(h,i 
hs 



/l,2 
/2,2 



fl,n 
f2,n 



\ 



\fn,X fn 



(ii) 



fn,n J 



The resulting reduced density matrices, in terms of the 
spinless fcrmion operators, takes the form [26] 

n 

Pn = IT (^™ C « C ™ + (1 _ Ki)c n cD (12) 

where A n is the n th eigenvalue of the correlation matrix 
F n and the values A n and 1 — A ra are the eigenvalues of 
the reduced density matrices. 



IV. RESULTS: GGM VS BIPARTITE 
ENTANGLEMENT 

We now demonstrate the efficacy of the generalized ge- 
ometric measure in detecting the quantum phase transi- 
tion in the three-body interaction Hamiltonian (Eq. (3)). 
We begin by considering the infinite chain. Later, we will 
also consider the situation for finite ones. 



A. Infinite spin chain 

The two-site reduced density matrix p m . m +i (after tak- 
ing the partial trace over all other sites except m and 
777 + 1) of the ground state, can be expressed as 



(u ON 

z y 

y z 

\0 u) 



(13) 



where the elements can be calculated from the correlation 
functions: 





+ (^m^m+l)) 








- /l,2^ ^2 _ - 








- (O'm^m+l)) 


i 


/l,2) 2 , 


4 + 


T ((^m^m+l) + ( a 


/l,2, 





(14) 



where the averages are carried out in the ground state. 

We can now calculate the bipartite entanglement of 
the two-party reduced state, obtained from the ground 
state, by using entanglement measures like concurrence. 
The concurrence of a two-qubit state p2 is defined as 



C = max{0, Ai — A2 — A3 — A4}, 



(15) 



where the A^'s are the eigenvalues, with Xi being the 
largest, of the Hcrmitian matrix R = \J ^fp^pi^fp^. Here, 
= ® C7 y)P2( a y ® <J y)^ where the asterisk indicates a 
complex conjugation in the computational (i.e. <r z ) basis. 
The behavior of concurrence of the two-particle reduced 
density matrix of the ground state with respect to a is 
depicted in Fig. 1. The derivative of concurrence with 
respect to a shows discontinuities (see inset of Fig. 1) 
at a — 1 and a = tt/2, although the phase transition in 
the model occurs at a = 1. Hence it is clear that in this 
three-body interaction Hamiltonian, bipartite entangle- 
ment cannot faithfully capture the quantum cooperative 
properties of the many-body system. More precisely, it 
signals to a QPT which is not present in the quantum 
system [20]. 

Cooperative phenomena like quantum phase transi- 
tions in a many-body system occurs due to the interplay 
among a large number of (if not, all the) particles of the 
system. This is despite the fact that the couplings of 
a given particle of the system are typically (although, 
not always) with a few "nearby" particles. The effect is 
arguably like in entanglement swapping, where two parti- 
cles can get entangled even if they have never interacted 
in the past, but have interacted with two other particles 
respectively who in turn have interacted in the past [27]. 
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FIG. 1. (Color online.) The failure of bipartite entangle- 
ment to indicate QPT in spin Hamiltonians with three-body 
interactions. The figure represents the variation of concur- 
rence (C) with respect to a. The concurrence (in ebits) is 
represented on the vertical axis and a (dimensionless) is rep- 
resented on the horizontal axis. The sharp change of C at 
a — 1 corresponds to a QPT. There is a second sharp change 
at a — it/2, which however does not correspond to any QPT. 
The inset figure depicts the variation of the first derivative 
(C) of concurrence with respect to a. 



Since a large number of particles arc involved in creat- 
ing the quantum phase transition, it is natural to expect 
that it will be signaled by a multisitc physical property 
of the system, like a genuine multipartite entanglement 
measure. 

Towards this aim, let us now consider the behavior of 
the GGM, which is a genuine multipartite entanglement 
measure, in the ground state of the three-body interac- 
tion Hamiltonian given in Eq. (3). As we have discussed 
in Sec. II, to evaluate the GGM, one has to calculate the 
Schmidt coefficients of different bipartitions. Before cal- 
culating GGM, for case of calculation, let us introduce 
here a multiparty entanglement measure, based on the 
GGM, for an TV-party pure quantum state \tPn), as 



G n {\ip N )) = 1 - max{A^. rest } 



(16) 



where {A^. rest } denotes the set containing the squares of 
the Schmidt coefficients of the state |t/>at) in a partition 
of the N constituent parties into a bipartition of n and 
N — n parties. This is still a multiparty entanglement 
measure (which we refer to as the "n-th order geometric 
measure"), and can be used to reveal information about 
the GGM of the TV-party state. The latter holds due 
to the fact that the minimum of G„'s corresponds to 
the GGM of the multiparty state of N parties. Since 
the interactions in the Hamiltonian are between nearest 
two spins or between nearest three spins, it is plausible 
that the most interesting physics will be seen in the G n 
obtained from nearest n spins, and the same is referred 
to as G„ below. With the help of the Eqs. (9), (10), and 
(11), we are able to find the analytical form for different 
G„'s in the bipartitions. In the case of bipartition of 2 



FIG. 2. (Color online.) Generalized geometric measure de- 
tects quantum phase transition. The figure represents the 
trend of the n-th order geometric measures (Gn), on the ver- 
tical axis, with respect to a on the horizontal axis. The inset 
depicts the same figure, but by zooming in near the phase 
transition. The G n 's increases monotonically after the tran- 
sition and reaches 1/2 for large a. All quantities are dimen- 
sionless. 



spins versus the remaining spins ("2:rest"), we have 
1 1 



G 2 = 



2 7T 
1 1 



a < 1, 



a>l. 



(17) 



2 7ra 

Similarly for n = 3 and n = 4, we respectively have 
1 \/2 



Go = 



1 s/2 

2 7ra 



a < 1, 
a > 1, 



(18) 



and 



Gi 




13 



3tt 

1 2 - V4 + 9a 4 
na 37ra 3 



a < 1, 
a>l. 



(19) 



It is clear from the expressions of G2 , G3 , and G4 that 
whatever bipartition we consider, the minimum value of 
the set {G2, G3, ...}, which correspond to the GGM of 
the ground state, is a constant for a < a c = 1, and for 
a > a c , it always increases as can be seen from Fig. 2. 
Note that, the G n 's have the order G2 > G3 > G4 > . . ., 
and they merge in the asymptotic limit a — > 00, to reach 
the value 1/2. 

In the regime a < a c = 1, the GGM can be expressed 
in the form 



G n = 1/2 



(20) 



where the A n is a constant (with respect to a) which only 
depends on the n. The constant A n slowly increases with 
the increase in n (see Table I). Hence it clearly converges 
to 1/2 as n — >• 00. Therefore for a < a c = 1, G n 
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2 


0.318309 
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0.450158 


4 


0.488664 
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0.497669 
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0.499544 
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0.499913 


n — ¥ oo 
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TABLE I. The behavior of A n with respect to n. 
A/J 




FIG. 3. (Color online.) Energy gaps for finite systems. We 
plot the energy gaps, A/J, for different finite-size spin chains 
against the system parameter a. Both axes represent dimen- 
sionless quantities. 

as n — > oo. So, for a < 1, genuine multipartite entan- 
glement is absent while a finite bipartite entanglement is 
present in the system. 

Let us now consider the case when a ^ a c = 1. In 
this regime, the GGM (G„) depends on a, and can be 
written in the form 

G n = l/2-0(l/a). (21) 

We find, from the expression of F n defined in Eq. (11), 
that as a — > oo, G n — > 1/2 irrespective of n. So, when 
a — > oo, the system becomes fully multiparty entangled 
with a vanishing bipartite entanglement. 

B. Finite spin chains 

Recent experimental developments strongly indicate 
that the Hamiltonians considered in this paper can be 
experimentally observed and the behavior of entangle- 
ment can be studied in a real system. For example, it 
has already been proposed that by using two-component 
ultracold Bose-Bose and Fermi- Fermi gas mixture and by 
suitably tuning the scattering length, the XX Hamilto- 
nian can be realized in ultracold atoms [3, 7]. Moreover 
the additional three-body interaction has also been cre- 
ated in the laboratories [19]. 

Keeping such possibilities in mind, in this section, we 
study the behavior of genuine multipartite entanglement 
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FIG. 4. (Color online.) The GGM for finite-size systems. 
We plot the GGM for systems of different finite spins (= N) 
on the vertical axis, against the system parameter a on the 
horizontal one. Both quantities are dimensionless. See text 
for the finite-size scaling analysis for the transitions. 

in controllable finite spin systems, which are achievable 
by using current technologies. The Hamiltonian that we 
consider here is same as in Eq. (3). We perform exact 
diagonalization to explore the properties of the ground 
state of the system. The spin chains with an odd number 
of spins have doubly degenerate ground states, while that 
is not the case for one with an even number of spins. The 
degeneracy in the case of an odd number of spins can be 
removed with the introduction of a small transverse field. 

The ground states have a vanishing total magnetiza- 
tion for the whole range of a. The energy difference (A) 
between the ground state and the first excited state has 
been calculated and the variation is shown in Fig. 3. 
In the thermodynamic limit, a quantum phase transition 
in general occurs with the vanishing of the energy gap. 
However, finite size simulations in our case show that 
there are several values of a where the A — > 0. It turns 
out that the derivatives of genuine multipartite entan- 
glement measure (GGM) also show discontinuities at the 
points where A — > 0, i.e. at points additional to the one 
near a = 1 (see Fig. 3). 

The derivative of the GGM exhibits a discontinuity 
near a = with approaching to 1 with an increase 
in the number of spins. This is the signature of the quan- 
tum phase transition at a — a c = 1 in the infinite chain. 
By performing a finite-size scaling, we find that the point 
of discontinuity present in the finite system, approaches 
to a c = 1 as 

a ~ a c + N- 1 - 787 . (22) 

Here N is the number of spins. There appears second dis- 
continuities (for different N), seen in Fig. 4, and these are 
away from a = 1 for relatively small N. For example, for 
N = 8, 10, and 12 the second ones occur at a = 2.6, 1.7, 
and 1.4 respectively It is clear that the second disconti- 
nuity point approaches to a = 1 in the large system limit. 
Again a finite-size scaling is employed, and we find that 
the point of the second discontinuities (a^ ) approaches 
to a c (= 1) as 

o% - a c + N-' 3A . (23) 
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V. CONCLUSION 

We have investigated the advantage of using the gen- 
uine multiparty entanglement measure called generalized 
geometric measure for detecting quantum phase transi- 
tions in the infinite quantum spin- 1/2 chains with two- 
spin isotropic XY interactions and three-spin interac- 
tions. We find that in contrast to bipartite entanglement, 
like concurrence, which signals to phase transitions that 
are not present, the generalized geometric measure faith- 
fully detects the quantum phase transition in the sys- 
tem. Recent experimental achievements that have made 
it possible to realize quantum spin models in a controlled 
way, motivated us to consider the behavior of the genuine 
multiparty entanglement for finite spin chains. We reveal 



that the generalized geometric measure can successfully 
notice changes in the finite systems, by analyzing the 
scaling behavior of which, we can detect the quantum 
phase transition in the systems. 
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